Conserving many body approach to the fully screened, infinite U Anderson model. 
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Using a Luttinger Ward scheme for interacting gauge particles, we present a conserving many 
body treatment of a family of fully screened infinite U Anderson models that has a smooth cross- 
over into the Fermi liquid state, with a finite scattering phase shift at zero temperature and a 
Wilson ratio greater than one. We illustrate our method, computing the temperature dependence 
of the thermodynamics, resistivity and electron dephasing rate and discuss its future application to 
non-equilibrium quantum dots and quantum critical mixed valent systems. 

PACS numbers: 73.23.-b, 72.15.Qm, 72.10.Fk 



The infinite U Anderson model is a central element in 
the theory of magnetic moments, in their diverse man- 
ifestations within antiferromagnets, heavy electron sys- 
tems and quantum dots. While the underlying physics 
of the Anderson impurity model is well understood, with 
a wide variety of theoretical techniques available for its 
description d], there are many new realms of physics that 
relate to it, such as quantum critical mixed valent and 
heavy fermion materials [3J or voltage biased quantum 
dots|J| where our tools and understanding of equilibrium 
impurity physics are inadequate. 

These considerations motivate us to seek new theo- 
retical tools with the flexibility to explore the physics 
of the Anderson model in a lattice and non-equilibrium 
environment. One of the well-established tools of many 
body physics is the Luttinger Ward (LW) scheme 5] . This 
scheme permits a systematic construction of many-body 
approximations which preserve important interrelation- 
ships between physical variables, such as the Luttinger 
sum rule, or the Korringa relationship between spin re- 
laxation and spin susceptibility, that are implied by con- 
servation laws. Recent work on the application of the LW 
scheme to gauge theories provides new insight into 
how such conserving approximations can be extended to 
encompass Fermi liquids in which ultra-strong interac- 
tions are replaced by constraints on the Hilbert space. 

In this paper we show how these insights can be applied 
to the infinite U Anderson model to obtain a description 
of the smooth cross-over from local moment to Fermi liq- 
uid behavior, in which the important conserving relation- 
ships between the physical properties of the ground-state, 
such as the Friedel sum rule and the Yamada-Yosida 
relationships between susceptibilities and specific heat, 
are preserved without approximation. Our method is 
scalable to a lattice, and can also be extended to non- 
equilibrium settings. To demonstrate the method we 
compute the temperature dependent thermodynamics, 
resistivity and electron dephasing rate in a single impu- 
rity model, and discuss how the method can be extended 
to encompass a non-equilibrium environment. 



Our starting point is a Schwinger boson representation 
for the SU(N) infinite U Anderson model. Schwinger 
bosons can be used to describe magnetically ordered local 
moments, and in the limit of large N, this description 
becomes exact 0- Recent progress [1,E3 has shown how 
this scheme can be expanded to encompass the Fermi 
liquid physics of the Kondo effect, through the study 
of a family of "multichannel" Kondo model, in which 
the number of channels K commensurate with rib = 25, 
(K = 25), where 5 is the spin of the local moment. The 
corresponding infinite U Anderson model is 
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Here cl creates a conduction electron with momentum 

kva 

k, channel index v G and spin index a £ 

The bilinear product b^Xv between a Schwinger boson 
and a slave fermion (or holon) operator \u creates a 
localized electron in the constrained Hilbert space, which 
hybridizes with the conduction electrons. The \N de- 
nominator in the hybridization ensures a well-defined 
large N limit. The energy of a singly occupied impu- 
rity eo is taken to be negative. The conserved operator 
Q = J2 a b\pa + XlXv replaces the no-double occu- 
pancy constraint of the infinite U Anderson model, where 
Q = 25 — K is required for perfect screening. 
The characteristic "Kondo scale" of our model is 
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where D is the conduction electron half bandwidth, T — 
npV 2 is the hybridization width and p is the conduction 
electron density of states. This relationship follows from 
general leading logarithmic scaling of the model. The 
non-trivial prefactor in the Kondo temperature depends 
on the third order terms in the beta function. 
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The next step is to construct the LW functional Y[G]. 
The variation of Y[G] with respect to the full Green's 
functions of the conduction, Schwinger boson and 
slave fermion fields, = SY/SG^, self-consistently de- 
termines the self-energies of these fields0,0- Y[G] is 
equal to the sum of all two-particle irreducible free energy 
Feynman diagrams, grouped in powers of 1/iV (Fig^a)). 
In our procedure we neglect all but the leading order 
O(N) term in the LW functional, and then impose the 
full self-consistency generated by this functional. 

In the formal large N limit, the conduction electron 
self-energy is of order 0(1/N), and at first sight, should 
be neglected. However, the effects of these terms on the 
free energy are enhanced by the number of channels K 
and spin components N, to give a leading order O(N) 
contribution to the free energy: it is precisely these terms 
that produce the Fermi liquid behavior|lflj. The com- 
plex issue of exactly when the conduction electron self- 
energies should be included is naturally resolved by using 
the leading order LW functional, subsequently treating 
1/N as a finite parameter, using the conduction electron 
self-energies inside the self-consistency relations. The ex- 
plicit expressions for the self-energies are then 
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where G b , G x and G c are the fully dressed imaginary 
time Green's functions of the boson, holon and conduc- 
tion electrons respectively. Equations J3J| are solved self 
consistently, adjusting the chemical potential A to satisfy 
the averaged constraint (Q) = K. 

A key step in the derivation of the Ward identities is 
the presence of a scale in the excitation spectrum, which 
manifests itself mathematically in our ability to replace 
Matsubara summations by integrals at low temperatures: 

Er duj 

This transformation is the key to derivations of the 
Friedel sum rule, the Yamada-Yosida and Shiba rela- 
tionships between susceptibilities, specific heat and spin 
relaxation rates [l^. Il3l Il4| . There is long history of 
attempts to develop conserving approximations to the 
infinite U Anderson model |l5|, [l6j, but each has been 
thwarted by the development of scale-invariant X-ray sin- 
gularities in the gauge or slave particle spectra for which 
the above replacement is illegal. In the multi-channel 
approach, the gauge particles - the formation of a stable 
Kondo singlet when 2S = K leads to the formation of a 
gap in the spinon and holon spectrum [see Fig. dc)], 
which permits a Fermi liquid spectrum to develop at 
lower energies. This feature permits all of the manip- 
ulations previously carried out on the finite U Anderson 
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FIG. 1: (a)LW functional grouped in powers of 1/iV. Solid 
line - Schwinger bosons, dashed lines - holons, double lines - 
conduction electrons. Loops contribute O(N), vertices 1/yN . 
(b) Circles: phase shifts extracted from the calculated t- 
matrices, for K = 1, n x m 0.04. Solid line stands for 
(K — n x )/(NK), dashed line shows the phase shift without 
imposing self consistency on the conduction electron propaga- 
tors, (c) The boson spectrum is gapped at low temperatures. 



model to now be extended to infinite U. Moreover, these 
relationships are all satisfied by the use of the leading 
LW functional. 

There are many important conserving relationships 
that are preserved, provided the spinons and holons de- 
velop a gap. The first of these is the Friedel sum rule (or 
in the lattice, the Luttinger sum rule), according to which 
the sum of the conduction electron phase shifts must be 



equal to the total charge K 



on the impurity, whith 



as the ground-state holon occupancy, so 
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Fig. 1 . contrasts the dependence of phase shift obtained 
from the t-matrix in our conserving scheme, with that 
obtained using E c from the leading order 1/N expansion. 
The relationship Q is obeyed at each value of N in our 
approximation. 

The quasiparticle spectrum in the Fermi liquid ground- 
state is intimately related to the conduction electron 
phase shifts via Nozieres Fermi liquid theory [J^- The 
change in the scattering phase shift in response to an ap- 
plied field, or a change in chemical potential, is directly 
related to the spin, charge and channel susceptibilities, 
which in turn, leads to a generalized "Yamada-Yosida 
identity" , distributing the NK degrees of freedom among 
the spin, flavor and charge sectors: 
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Here 7 stands for the specific heat coefficient 7 — C/T, 
and %sj Xf an( l Xc stand for the spin, flavor and charge 
susceptibilities. Xs,/,c denotes the corresponding suscep- 
tibilities in the absence of the impurity. Our identity (0 
reduces to the known results for the Wilson ratio in the 
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FIG. 2: Upper pannel: The impurity contribution to the 
entropy vs. temperature. Inset: at low temperatures the en- 
tropy becomes linear and 7 saturates. Lower pannel: impurity 
susceptibility vs. temperature. For parameters see Fial3l 




FIG. 3: t-matrix for eo/D = -0.2783, T/D = 0.16, N = 4 
and K = 1 (T K /D = 0.002). Main frame: T/D=0.1, 0.08, 
0.06, 0.04, 0.02, and 0.01. Inset: T/10~ 4 D=10, 8, 6, 4, 2, 1, 
and 0.5. The dot marks sin 2 {7r(i\ - n x )/NK} 



one channel model^j| and the Kondo limit^^ ■ Since our 
scheme preserves the Friedel sum rule for each channel, 
it automatically satisfies this relationship. 

Various thermodynamic quantities can be extracted 
from the numerical solution of the self-consistency equa- 
tions. The development of a fully quenched Fermi liquid 
state is succinctly demonstrated by calculating the en- 
tropy, using a formula of Coleman, Paul and RechQ 
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where the trace denotes a sum over the spin and chan- 
nel indices of the various fields and tlq is the distribution 
function (Fermi for £ = c, x and Bose for £ = b). Gq and 
respectively denote the retarded Green functions and 
self energies. The single and double primes mark the real 
and imaginary part respectively. The impurity contribu- 
tion to the entropy <Si mp = S—Sq, where So is the entropy 
of the bare conduction band, is shown in the upper pan- 
nel of Fig. [21 An inspection of low temperatures shows 
a linear behavior of the impurity contribution to the en- 
tropy. The latter marks the saturation of the specific 
heat coefficient 7i mp = dS lmp /dT at low temperatures 
and the formation of a local Fermi liquid. 

To calculate the spin susceptibility, we couple a small 
magnetic field to the magnetization, defined as 
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The impurity magnetization is obtained by subtracting 
the magnetization of the free conduction sea from M. 
When we apply a field, we must keep track of the ver- 
tex corrections that arise from the field-dependence of 
the self-energies. These vertex functions are essential to 
maintain the conservation laws. We compute the mag- 
netic field vertex functions by iterating the solutions to 



self-consistency in a small magnetic field 11]. Figure © 
shows the impurity susceptibility x{T)- The susceptibil- 
ity is peaked slightly below the Kondo temperature and 
saturates to a constant at a lower temperature. 

One of the important conservation laws associated 
with the conservation of spin, is the Korringa Shiba re- 
lation between the dynamical and static spin suscepti- 
bility. On the assumption that the gauge particles are 
gapped, the LW derivation carried out by Shiba on the 
finite U Anderson model, some thirty years ago, can 
be simply generalized to the infinite U model, to ob- 
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guarantees that the power-spectrum of the magnetiza- 
tion is linear at low frequencies, ensuring that the spin 
response function decays as 1/t 2 in time. Remarkably, a 
Fermi liquid, with slow gapless spin excitations is sand- 
wiched beneath the spinon-holon continuum. 

We now turn to a discussion of the electron scatter- 
ing off the impurity, which is determined by the conduc- 
tion electron t-matrix, t(z) = - — ~(Sj^T: where G° c (z) = 

^2^{z — e c) 1 is the bare local conduction electron prop- 
agator. Various important physical properties, such as 
the impurity resistivity and the electron dephasing rate 
in the dilute limit can be related to this quantity. Figure 
13 shows the imaginary part of t(e + iff): as temperature 
is lowered towards the Kondo temperature, a resonance 
peak starts to develop around the Fermi energy. The 
peak continues to evolve with temperature and finally 
saturates a decade below the Kondo temperature. The 
fully developed resonance is pinned to the Langreth sum- 
rule [13 value of sin 6 c /ttp, at the Fermi energy presented 
by the dot in the inset of Fig02 

The impurity contribution to the resistivity R t is re- 
lated to a thermal average of the t— matrix, 
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FIG. 4: The impurity contribution to the resistivity Ri and 
dephasing rate t^ 1 as a function of temperature. Solid lines: 
the parameters used in Fig. |3] dashed lines: F = 0.1D. 



where hp is the Fermi wavelength and / is the Fermi 
distribution function. Figure 01 shows Ri(T). The resis- 
tivity increases as the temperature decreases, saturating 
at a value determined by the scattering phase shift. 

As a final application of the t-matrix, we compute the 
electron dephasing rate t^ 1 that controls the field depen- 
dence of weak electron localization. In the cross-over to 
the Fermi liquid state, the formation of the Kondo reso- 
nance gives rise to a peak in the inelastic scattering 
and the electron dephasing rate 20]. In the dilute limit, 
the impurity contribution to the dephasing rate is given 
by r^ 1 = 2ni[-Imf - np\t\ 2 } [H El H3] . Fi S- El shows 
the dephasing rate on the Fermi surface computed from 
t, showing how the scattering scales with the Kondo tem- 
perature Tff in the Kondo regime of the model. 

Our calculation of the dephasing rate shows that our 
method captures the cross-over into the coherent Fermi 
liquid. However, it also highlights a shortcoming that we 
hope to address in future work. At low temperatures, the 
electron dephasing rate has a quadratic dependence on 
temperature. The electron-electron scattering diagrams 
responsible for these processes involve an internal loop of 
gauge particles, which only enters in the subleading O(l) 
terms of the LW functional shown in Fig. 1(a), which 
are absent from the current work. On the other hand, 
the low temperature, on-shell value of the conduction 
electron vertex is directly related to 7 by Ward identities 
[2ll | , providing a possible future simplification for treating 
these processes. 

There are many directions for future development. 
We are particularly interested in the extension to non- 
equilibrium quantum dotsQ. The theoretical under- 
standing of the Kondo physics of quantum dots at finite 
voltage bias is still evolving. An extension of our method 
to a voltage biased infinite Anderson model can be made 
using a Keldysh generalization of the self-consistency 
equations |23T|. Many features of our approach, its con- 
serving properties, its inclusion of scaling properties up 
to third order in the beta function, and its systematic in- 



corporation of spin relaxation effects on the boson lines, 
suggest that this will be a robust method to examine how 
spin and electron dephasing effects evolve with voltage in 
quantum dots. These methods can also be used to study 
the voltage dependence [24| of the distribution function 
and noise in magnetically doped mesoscopic wires |25|. 

One of the striking features of the Schwinger boson 
approach to the Anderson and Kondo models, is the co- 
existence of a gapless Fermi liquid, sandwiched at low 
energies between a gapped spinon and holon fluid. The 
gap in the spin-charge decoupled excitations appears inti- 
mately linked to the development of the Fermi liqid Ward 
identities. The methods presented in this paper can be 
scaled up to describe the dense Anderson and Kondo lat- 
tice models, where the Friedel sum rule is replaced by 
the Luttinger sum ruleQ. In the lattice, the U(l) lo- 
cal symmetry of the impurity model will in general be 
replaced by a global U(l) symmetry associated with the 
pair-condensed Schwinger bosons 26J. In this setting, the 
gapped spinons and holons are propagating excitations, 
whose gap is fundamental to the large Fermi surface. It 
is this very gap that we expect to collapse at a quan- 
tum critical point. The current work, which includes the 
effects of valence fluctuations, provides a powerful frame- 
work for examining this new physics. 
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